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Abstract 

Nonlinear waves are studied in a mixture of liquid and gas bubbles. 
Influence of viscosity and heat transfer is taken into consideration 
on propagation of the pressure waves. Nonlinear evolution equations 
of the second and the third order for describing nonlinear waves in 
gas-liquid mixtures are derived. Exact solutions of these nonlinear 
evolution equations are found. Properties of nonlinear waves in a 
liquid with gas bubbles are discussed. 



1 Introduction 

A mixture of liquid and gas bubbles of the same size may be considered as an 
example of a classic nonlinear medium. In practice analysis of propagation of 
the pressure waves in a liquid with gas bubbles is important problem. Similar 
two - phase medium describes many processes in nature and engineering 
applications. In particular, such mathematical models are useful for studying 
dynamics of contrast agents in the blood flow at ultrasonic researches [HE]- 
The literature on this subject deals with theoretical and experimental studies 
of the various aspects for propagation of the pressure waves in bubbly liquids. 

The flrst analysis of a problem bubble dynamics was made by Rayleigh [3j , 
who had solved the problem of the collapse of an empty cavity in a large mass 
of liquid. He also considered the problem of a gas - fllled cavity under the 
assumption that the gas undergoes isothermal compression [1]. Based on 
the works of Rayleigh [3] and Foldy [5] van Wijngarden [B] showed that in 
the case of one spatial dimension, the propagation of linear acoustic waves 
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in isothermal bubbly liquids, wherein the bubbles are of uniform radius, is 
described by the linear partial differential equation of the fourth order [7]. 

The dynamic propagation of acoustic waves in a half - space filled with 
a viscous, bubbly liquid under van Wijngaarden linear theory was consid- 
ered in the recent work [7j. However we know that there are solitary and 
periodic waves in a mixture of a liquid and gas bubbles and these waves can 
be described by nonlinear partial differential equations. As for examples of 
nonlinear differential equations to describe the pressure waves in bubbly liq- 
uids we can point out the Burgers equation [814T0]. the Korteweg - de Vries 
equation [HHIl], the Burgers - Korteweg - de Vries equation [H] and so on. 

Many authors applied the numerical methods to study properties of the 
nonlinear pressure waves in a mixture of liquid and gas bubbles. NigmatuUin 
and Khabeev [12] studied the heat transfer between a gas bubble and a 
liquid by means of the numerical approach. Later Aidagulov and et al. |16j 
investigated the structure of shock waves in a liquid with gas bubbles with 
consideration for the heat transfer between gas and liquid. Oganyan in [T7p8] 
tried to take into account the heat transfer between a gas bubble and a liquid 
to obtain the nonlinear evolution equations for the description of the pressure 
waves in a gas-liquid mixture. However, we have some different characteristic 
times of nonlinear waves in these processes and it turned out that the solution 
of this task is difficult problem. 

The purpose of this work is to obtain more common nonlinear partial 
differential equations for describing the pressure waves in a mixture liquid 
and gas bubbles taking into consideration the viscosity of liquid and the 
heat transfer. We also look for exact solutions of these nonlinear differential 
equations to study the properties of nonlinear waves in a liquid with gas 
bubbles. 

This Letter is organized as follows. System of equations for description 
of nonlinear waves in a mixture of liquid and gas bubbles taking into con- 
sideration for the heat transfer and the viscosity of liquid is given in section 
2. In section 3 we obtain the basic nonlinear evolution equation to describe 
the pressure waves in liquid with gas bubbles. In sections 4 and 5 we present 
nonlinear evolution equations of the second and the third order and search 
for exact solutions of these nonlinear differential equations. 
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2 System of equations for description of mo- 
tion of liquid with gas bubbles with consid- 
eration for heat exchange and viscosity 

Suppose that a mixture of a liquid and gas bubbles is homogeneous medium 
[T9] . In this case for description of this mixture we use the averaged temper- 
ature, velocity, density and pressure. We also assume that the gas bubbles 
has the same size and the amount of bubbles in the mass unit is constant 
A^. We take the processes of the heat transfer and viscosity on the bound- 
ary of bubble and hquid into account. We do not consider the processes of 
formation, destruction, and conglutination for the gas bubbles. 

We have that the volume and the mass of gas in the unit of the mass 
mixture can be written as 

V=^nR^N, X = Vpg, 

where R = R{x,t) is bubble radius, Pg = Pg{x,t) is the gas density. Here 
and later we believe that the subscript g corresponds to the gas phase and 
subscript / corresponds to the hquid phase. 

We consider the long wavelength perturbations in a mixture of the liquid 
and the gas bubbles assuming that characteristic length of waves of pertur- 
bation more than distance between bubbles. We also assume, that distance 
between bubbles much more than the averaged radius of a bubble. 

We describe dynamics of a bubble using the Rayleigh - Lamb equation. 
We also take the equation of energy for a bubble and the state equation for 
the gas in a bubble into account. The system of equation for the description 
of the gas bubble takes the form P^l^ 

ft (flfl,, + + ^fl,] = -p, - -p, (1) 



r. = ^f#|, (3) 



Pg,0 \Ro 

where P{x,t) is a pressure of a gas-liquid mixture, Pg is a gas pressure in 
a bubble, Tg and Ti are temperatures of liquid and gas accordingly, Xg is a 
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coefficient of the gas thermal conduction, Nu is the Nusselt number, n is a 
polytropic exponent, v is the viscosity of a hquid. 

The expression for the density of a mixture can be presented in the form 

m 

.^—,V^,^—±—. ,4, 

Considering the small deviation of the bubble radius in comparison with 
the averaged radius of bubble, we have 

R{x,t) = Rq + rj{x,t), Rq = const, ||?7||<<_Ro! R{x,0) = Rq. (5) 

Assume that the liquid temperature is constant and equal to the initial 
value 

Ti = T \t=o= To, To = const. 
At the initial moment, we also have 

t = 0: P = Pg = P^^ Pq = const, V = Vq = -nR^N. 

o 

Substituting Pg and Tg from Eqs.([l]) and ([3]) into Eq.(j2]) and taking rela- 
tion (|5]) into account we have the pressure dependence of a mixture on the 
radius perturbation in the form 



W + + * qP ^" ^"r2 



, pii8uH{3n-l) + 9Rl) , Aupi 2P, , 



2RlP, 



?,xNu{n - l)To' 

(6) 

From Eq.dl]) we also have the dependence p on rj using formula (|5]) 

2 P« 
P = PO - pi] + piT] , po 



1-X + Vopi 

(7) 

Sp^Vo _ QpfVo{2piVo-l + X) 

^ Roil-X + Vopiy' Rlil-X + piVor ■ 



We use the system of equations for description of the motion of a gas- 
hquid mixture flow in the form 



dt dx ' ^ \ dt ^ dx ) dx 

where u = u{x, t) is a velocity of a flow of a gas-liquid mixture. 

Eq.® together with Eqs.Q and ([7]) can be applied for description of 
nonlinear waves in a gas-liquid medium. 

Consider the linear case of the system of equations ([7]) and (jS])- As- 
suming, that pressure in a mixture is proportional to perturbation radius, 
we obtain the linear wave equation for the radius perturbations 



Vtt = clr]^^, co = y-^. (9) 
Let us introduce the following dimensionless variables 

t = — t', X = lx', U = Cou', T] = RoT] , P = PoP' + Po, 

Co 

where I is the characteristic length of wave. 

Using the dimensionless variables the system of equations ([6]), ([7]) and ([8]) 
can be reduced to the following (the primes of the variables are omitted) 

fiRo p 

^° {ut + uu^) + riut - \Px = 



/ii?o ' ' ' 3 (10) 

P + x,Pt + r]P + 3nx, r]tP = -(A + /^s) Vtt - {W2 - A) mtt- 



2 2^ 

where the parameters are determined by formulae 

xpiRlcl XCo 
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3 Basic nonlinear evolution equation for de- 
scription of nonlinear waves in liquid with 
bubbles 

We can not find the exact solutions of tlie system of nonlinear differential 
equations (|T0|) . To account for the slow variation of the waveform, we intro- 
duce a scale transformation of independent variables [21] 

^ = e"'{x-t), T = e'^+H, m>0, e < 1, (12) 

dx d^' dt dr 

Substituting (IT^ into (ITU]) and dividing on e"^ in first two equations we 
have the following system of equations 

^Vt -Vi ^ui. + rjUf: + uri(. - e rjr]^ H rjr]^ = 0, (13) 

f Po \ , Po , Po Ion 

' ['Wo ) + Wo"' + '""^ ~ ~ " 3 = °' ^'^^ 

P + e^+^xiP^ - e'^xiP^ + r/P + e^+^Snxir/^P - e^SnxiT^^P = 

(2/32 - /3i) Wr. + £'™+'2 (2/32 - /3i) - 
- (2/32 - /3i) - e^'"-^^ + ^Z^^) r/^+ (15) 

+^^-"^2 (^/3, + ^/32) - (^/3, + ^/32) ,|- 

We now assume that the state variables u, r/ and P can be represented 
asymptotically as series in powers of e about an equilibrium state 

M = £Mi + E^U2 + . . . , r] = er]i + e^r/a + • • • , P = ePi + e^Pa + . . . (16) 
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Substituting (fT^ into flT^- ffT^ and equating expressions at e and e"^ to 
zero, we obtain some equations with respect to Ui, M2, ?7i, rj2, Pi and P2. 
Solving these equations with respect to Pi we have the equation for pressure 
Pi in the form 

^ir + a nn^ + £ — ^ — n^se - e — — — PiP^h - 



where 



a = + -. 

Po p 3 



From Eq. f ITTj) one can find some nonhnear evolution equations for de- 
scription of waves in a mixture liquid and gas bubbles with consideration for 
the heat transfer and the viscosity. 



4 Nonlinear evolution equation of second or- 
der for description of waves in liquid with 
bubbles 

Assuming m = 1, (/3i + /32)/6 = /^e'', 5 > 0, /3 ~ 1 in Eq.([I7]), we have 
nonlinear evolution equation of the second order in the form 

Pi. + a PiPi, = - f) Pi,i + e ^ {PiPid, . (19) 

Using transformations Pi = 1/e P[ and a = ea one can write Eq. f[TI?|) in 
the form 

P[^ + aP[p[^ = AP[^^ + Ai {p[Pu) , (20) 



A _ xi ^ Co / _ PqRI \ 
6 2 ~ 3/ V3Po xNuTo) ' 

nxi n Co Pq 

~2~ ~ 3xNu{n-l)Tol' 



We can see that dissipation of the nonhnear waves is described by the 
equation f l2U]) and depends on values of coefficients A and Ai. 

The first term in expression for A corresponds to a dissipation of a wave 
because of the viscosity on boundary a bubble-hquid. The second term in 
A corresponds to the heat transfer between a hquid and bubbles. We have 
to note that there is the decrease of the coefficient A in the case of the heat 
transfer. The coefficient Ai characterizes the nonlinear dissipation of a wave 
which is explained by the heat transfer between liquid and gas bubbles. 

In the case of the isothermal process {Nu — )■ cxd, xi — )■ 0) Eq. fl20|) goes to 
the Burgers equation. We have the dissipation of the wave that is determined 
by the value of the viscosity of a liquid. Currently the Burgers equation is 
well known. The wave processes described by this equation are well known. 
Taking the Cole - Hopf transformation into account this equation can be 
reduced to the heat equation [11, 12]. Analytical solutions of many problems 
described by this equation were constructed. 

The nonlinear term in a right hand side (1201) leads to the additional damp- 
ing of the pressure wave in comparison with the Burgers equation. These 
nonlinear dissipative effects are connected with the heat transfer in a mix- 
ture liquid and gas bubbles. 

One can construct some exact solutions of the equation fl20l) . Using the 
scale transformations in Eq. (l20|l 



^ 2a' ^ ' (A-3xi) V 2a' / ' ^ 3nxi ' ^ ^ 

we obtain the following equation (the primes are omitted) 

Vr +VVi:= Vj:^ + {V Vj:)^ . (22) 

Some exact solutions of Eq. (l22|) can be found. Taking transformation 
(pT]) into account we obtain the solution of Eq. (!20|) by the formula 

-1 

± , " I -1 



3n xi 



n Xi 6 \n Ki J 

(23) 

From expression ( 125]) one can see, that the amplitude of a wave depends 
on its velocity. Besides the amplitude of this pressure wave is defined by a 
relation between viscosity of a liquid and a gas thermal conduction. 



8 



'(1) 

Dependencies of solution from ^ at different values parameter xi 
at Co = 2, C2 = 1, r = are illustrated on Fig.l. We use values xi are 
following: 0.200; 0.205; 0.210. From Fig.l we can see that the amplitude of 
a wave decreases with growth xi. Similar profiles of the pressure waves were 
observed experimentally. They are called by shock waves with monotone 
structure or waves 



Pi 1 




Figure 1: Kink wave solution 022]) of Eq. ([20]) at x = 0, 20; 0, 22; 0, 24 
(curves 1, 2, 3) 



Experimental study of the damping of the pressure waves in a liquid 
with a mixture of bubbles of two different gases were presented in |25]. In 
particular it was found the damping of amplitude for the pressure wave for 
a liquid with freon bubbles. The volume gas content of freon bubbles was 
a = 1.2% and the diameter of bubbles was d = 2,3 mm. 

We have considered the propagation of the damping of the pressure im- 
pulse which is described by Eq.( l20|) . For numerical modeling of the solitary 
waves we use the difference equation which corresponds to Eq. (l20|) . We have 
observed the good agreement of results of numerical modeling with experi- 
mental data. This comparison is given on Fig. 2. 

Let us discuss the features of the wave dynamics described by the equation 
fl20|l . For numerical modeling of the solitary waves we use the difference 
equation which corresponds to Eq. fl20|) . 

Comparison of the solitary pressure pulses described by the Eq. fl20l) with 
the solitary pressure pulses of the the Burgers equation ((Ai = in fl2Ul) )is 
given on the Fig. 3. We can see from Fig. 3 that pressure front described by 
Eq. fl20l) is more smooth than pressure front described by Burgers equation. 
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Figure 2: Comparison of amplitude attenuation of a pressure wave described 
by the equation fl20|l with experimental data 



So the heat transfer between the gas in bubbles and the liquid leads to 
additional dissipation of the solitary wave. 




Figure 3: Evolution of the pressure solitary wave described by Burgers equa- 
tion (left figure) and Eq. (I2U]) (right figure). 
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5 Nonlinear evolution equation of the third 
order for description of waves in a liquid 
with gas bubbles 

Assuming that m = i, A = Ge^'^/^A', A' ~ 1, xi = 26^/^ x^, ~ 1 in 

Eq. (1171) and using transformations Pi = — P[, a = e a we have the nonhnear 
evolution equation of the third order 

P' +aP'p' +l^l±Ap' 2/32 - 1^1 p'p' _ 



The coefficient of dispersion can be determined taking the relations f jTTj) 
into account 

6 6/2 Po V 9x,iVn(n-l)Toy'- ^ ^ 

We can see that the value of dispersion depends on the initial radius of 
a bubble and a relation between a viscosity and the coefficient of the heat 
transfer on interphase boundary. In the case of the isothermal condition for 
gas in bubbles {Nu — )■ oo) or in the case of small viscosity of a liquid (z/ — )■ 0) 
the wave dispersion is defined by value of radius of bubbles. 

Taking Nu — ?■ oo and z/ — )■ into account in Eq. fl2^ we have Korteweg 
- de Vries equation that was obtained for description of nonlinear waves in 
a liquid with gas bubbles [H]. Using Nu — )• oo we have from Eq.( l2^ the 
Burgers - Korteweg - de Vries equation, obtained in |8]. This equation is 
not integrable but some special solutions of Burgers - Korteweg - de Vries 
equation were obtained in |221En]- It is likely that the nonlinear evolution 
equation flMj) is new for the description of nonlinear waves in a liquid with 
gas bubbles. 

Consider Eq. flMj) without a dissipation. In this case we have 
A + /32 2/32 - /3i 

^ (26) 



3n-2 ^ 5 



18 - ■ 18'-/ 
Using the scale transformations 

r-^—r' p'-^-i^l±M, (27) 
^"/3i + /3/' - 2/32-/31 ^' ^'^^ 
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we obtain equation (the primes are omitted) 

9a' Sla'PQl'^xNu{n-l)TQ 



2/32 - /3i 2pi cl Rli'^xNuin - l)To - 4z/Po) ' 



(3n - l)/3i + 3/32 _ 27x Nu{n - l)To + 8 Nu{?,n - 1)Pq 

^ " 2/32 - /3i ~ 9xAr?x(n-l)To-4z/Po ' 

The travehng wave reduction of Eq.f l28p coincides with the Camassa - 
Holm equation |2I1E2] that has one of interesting class of solutions in the 
form of peaked solitons . The most simple form of these solutions is described 
by the formula 

v{^,t) = Aexp{-\k^-k^T + B\), (29) 
where k is wave number, A and B are arbitrary constants. Taking 

/5 = t-2, (30) 
and fg'g'g' = —v^'^'r' into consideration Eq. fl25]) can be written as 

nir' + V m^i + — 1 j v^f m = 0, (31) 



1 r+oo 

v{0 = 7^ / e-\^'-y\ m{y) dy, m = v- v^,^,, r' = k^ r, ^' = k^. (32) 

Using relations ( E7|) , ( 15U]) , we have the solution of Eq. (12^ in the form of 
the peaked solutions 

P^(,,.-,.m±^A^^,(.,,,JA±MJ^^,S\). (33, 

2/32 - Pi 6 

The graph of a solution at A = 1/6, S = 5, /32 = /3i = l,r = is 
demonstrated on Fig. 4. 

Using the traveling wave ansatz f(^,r) = y{z),z = ^ — Cqt in Eq. fl28|) 
and integrating on z, we obtain the equation 

Cl - Coy + + -yy.. - ^yj = 0. (34) 



12 



2 



6 



10 



5 14 



Figure 4: Peaked solution ([33]) of Eq. ([26]) 



The general solution of Eq. ([M]) at Co = ai, C\ = a\j2 is expressed by 
formulae 



y{z) = 1- fCse^^^^)'^" 



(35) 



where C3 7^ is constant of integration. Accordingly to fl35]) there are two 
families of solutions of Eq. ([3i 



VI«i(/3 + 2)| 



1 -l/(/3+2) 



ai(/3+2) < 0, 
(36) 



y(z) = 1- 



tanh 



V«i(/3 + 2) 



-l/(/3+2) 



;^+c^2) -1 



ai(/3+2) > 0, 
(37) 



where Ci is arbitrary constant. 

Family of solutions ai(/3 + 2) < describes the periodic waves. In the 
case /3 < — 2 solution (l36l) has the singular points that are determined by 
formula 



— —C2 + 



VI«i(/3 + 2)| 
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Family solutions at a;i(/3 + 2) > and /3 < —2 are solitary waves. In the 
case P > —2 solution goes to indefinite when z — )■ ±oo. Dependence of this 
solutions from z at parameter /3 = — 1 at C2 = 0, C3 = 1 is presented on Fig. 
5. 




Figure 5: Periodical solution (136|) of Eq. at /3 = — 1 



One can find some exact solutions of the equation (l24l) with the dissipa- 
tion. Using the scale transformations 

, 2(32 - (3i y 6 f 2(32-Pi Y ' p'_^(A+M c,^^ 

we obtain from equation (1241) the following equation (the primes are omitted) 

Vr + aivv^ + f = {v Vj:^)^ + iSv^v^i: + a Vj:^ + {y f 5)^ , (39) 

where 

^ (2/32 - ^ (3n - l)/3i + 3/32 ^ (V - x;)(2/32 - /3i) 

ISn^x;^ ' ^ 2/32 -/3i ' 3nx;(/3i + /32) ' 

Eq. fl5^ does not pass the Painleve test and is not integrable by the inverse 
scattering transform but there are some exact solutions of this equation. 
These solutions can be found using the simplest equation method [33ti36] for 
the travelling wave f (^, r) = y{z), z = ^ ~ Cqt. We look for solution in the 
form 

y{z) = ao + aiw{z) + a2w{z)'^, (40) 
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where w{z) is solution of linear ordinary differential equation of the second 
order [571155] 

w,, = Aw, + Bw + C. (41) 
One of these solution at ai < — 1 takes the form 
/ N , «i(«o -l)F f ai 

(42) 

ai 

where ag, C, C2 and C3 are arbitrary constants. 

Exact solution f H2]) of Eq. fl5I?]) represents oscillating fast damping waves 
which can be observed at propagation of nonlinear waves in a mixture liquid 
and gas bubbles. Analyzing the wave dynamics described by the nonlinear 
ordinary differential equation of the third order ( 139|) we can see that this 
equation describes the different forms of the solitary and periodical waves. 



63 sm < \ + C2 cos < M e 2 



6 Conclusion 

We have considered nonlinear waves in a mixture of a liquid and gas bubbles 
taking into consideration the viscosity of liquid and the heat transfer between 
liquid and gas bubbles. We have applied different scales of time and coor- 
dinate. As a result we have obtained the nonlinear evolution equations ( |20|) 
and ( 12^ of the second and the third order for describing the pressure waves 
in a liquid with bubbles. It is likely that these differential equations are new 
and generalize the Burgers equation, the Korteweg - de Vries equation and 
the Burgers - Korteweg - de Vries equation. These evolution equations allow 
us to take into account the influence of the viscosity and the heat transfer 
on the boundary of liquid and bubbles. Generally these equations are not 
integrable but all these equations have some exact solutions that have been 
found in this Letter. 

We had observed that nonlinear evolution equation f l2B]) can be used 
for describing waves at characteristic time t ~ Without dissipative 

processes Eq. (l26|) has a solution in the form of peaked solitons (so-called 
peakons). Taking into account dissipative processes we have solution of 
Eq. fl24p for description of nonlinear waves. Exact solutions have been found 
and we obtain the damping of the periodic pressure waves. To describe the 
pressure waves in the case of t ~ the nonlinear evolution equation of the 
second order can be used. Using the travelling wave ansatz the analytical 
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solutions of this equation have been obtained. These solutions correspond to 
shock waves with monotone structure. 
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